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ABSTRACT 

We present new theoretical estimates of the relative contributions of unresolved blazars and star-forming galaxies to 
the extragalactic y- ray background (EGB) and discuss constraints on the contributions from alternative mechanisms 
such as dark matter annihilation and truly diffuse y- ray production. We find that the Fermi source count data do 
not rule out a scenario in which the EGB is dominated by emission from unresolved blazars, though unresolved 
star-forming galaxies may also contribute significantly to the background, within order-of-magnitude uncertainties. 
Ir. addition, we find that the spectrum of the unresolved star-forming galaxy contribution cannot explain the EGB 
spectrum found by EGRET at energies between 50 and 200 MeV, whereas the spectrum of unresolved flat spectrum 
radio quasars, when accounting for the energy-dependent effects of source confusion, could be consistent with 
the combined spectrum of the low-energy EGRET EGB measurements and the Fermi- Large Area Telescope EGB 
measurements. 

Key words: galaxies: active - galaxies: general - gamma rays: diffuse background - gamma rays: galaxies 
Online-only material: color figures 


1. INTRODUCTION 

Studies of the extragalactic y-ray background (EGB) can 
provide insight into high-energy processes in the universe and, 
as such, has been the subject of much debate, particularly 
concerning the roles of extragalactic astrophysical sources and 
new physics . Recent data from the Large Area Telescope (LAT) 4 
on board die Fermi Gamma Ray Space Telescope allow for a 
reassessment of the possible astrophysical origins of the EGB, 
which could improve our understanding of y- ray production in 
these objects and provide more robust constraints on the more 
exotic scenarios. However, in order to determine the strength 
and spectrum of this isotropic background one needs to proceed 
from the raw photon count data by determining to the best extent 
possible the detector sensitivity, the intrinsic events produced by 
the larger charged particle flux impinging on the detector, and 
the much larger y-ray foreground within our Galaxy resulting 
from cosmic-ray interactions with photons and gas nuclei. 
Such analyses have been made for both the Energetic y-ray 
Experiment Telescope (EGRET) aboard the Compton Gamma 
Ray Observatory (Sreekumar et al. 1998; Strong et al. 2004a) 
and Fermi (Abdo et al. 2010j). 

Various extragalactic y-ray production scenarios have been 
explored theoretically as candidate components that could 
contribute significantly to the observed background. Among 
those considered are unresolved astronomical sources, such as 
active galactic nuclei (AGNs; Padovani et al. 1993; Stecker et al. 
1993; Salamon & Stecker 1994; Chiang et al. 1995; Stecker & 
Salamon 1996; Kazanas & Perlman 1997; Chiang & Mukherjee 
1998; Mukherjee & Chiang 1999; Miicke & Pohl 2000; Giommi 
et al. 2006: Narumoto & Totani 2006; Dermer 2007; Pavlidou & 
Venters 2008; Inoue & Totani 2009; Venters et al. 2009; Venters 
2010; Abazajian et al. 2010), star-forming galaxies (Pavlidou 
& Fields 2002; Fields et al. 2010; Makiya et al. 2011), and 
starburst galaxies (Thompson et al. 2007 ; Stecker 2007 ; Makiya 
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et al. 2011). The large majority of associated extragalactic 
sources detected thus far by both EGRET and Fermi are blazars 
(Hartmann et al. 1999; Abdo et al. 2010f), i.e., those AGNs 
for which the jet is closely aligned with the observer’s line 
of sight (Blandford & Konigl 1979), including y-ray loud flat 
spectrum radio quasars (FSRQs) and BL Lacertae type objects. 
It is expected that since blazars comprise the largest class of 
identified extragalactic y-ray sources, unresolved blazars should 
contribute significantly to the EGB. Additionally, just as our 
Galaxy produces y-rays, it is expected that y-rays are produced 
in other galaxies, and as such, unresolved galaxies might also 
contribute to the EGB with the most significant contribution 
originating from the population of actively star-forming galaxies 
(Stecker 1975; Pavlidou & Fields 2001, 2002; Fields et al. 2010; 
Makiya et al. 2011). Interesting truly diffuse mechanisms that 
could contribute to the EGB involve cosmic-ray interactions 
with intergalactic gas and the cosmic background radiation 
(Fazio et al. 1966; Stecker 1973; Dar 2007; Keshet et al. 
2003) and electromagnetic cascades produced by interactions of 
very high and ultrahigh energy particles with the extragalactic 
background light (Kalashev et al. 2009; Berezinsky et al. 2011; 
Ahlers et al. 2010; Venters 2010), as well as more exotic 
scenarios such as dark matter annihilation (Silk & Srednicki 
1984; Stecker et al. 1985; Rudaz & Stecker 1988; Stecker & 
Tyika 1989a, 1989b; Rudaz & Stecker 1991; Ullio et al. 2002) 
and decay (Olive & Silk 1985; Stecker 1986; Ibarra & Tran 
2008). 3 

In this paper, we estimate the contributions to the EGB 
from unresolved extragalactic y-ray sources of various types 
and compare them with the EGB obtained from analysis of 
Fermi data. In doing so, we also take into consideration the 
effects of both the completeness of the Fermi flux-limited blazar 
survey and the important effect of source confusion owing to the 
energy-dependent angular resolution of the Fermi-LAT detector. 
We will then briefly discuss the implications of possible truly 
diffuse emission mechanisms to the EGB. 

5 For reviews on dark matter annihilation, see Jungman et al. (1996) and 
Bertone et al. (2005). 
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Figure 1. Source count distributions for Fermi (black triangles; Abdo et al. 
2010f) and EGRET (red circles; Reimer & Thompson 2001) blazars. 

(A color version of this figure is available in the online journal.) 


2. EGRET AND FERMI RESOLVED SOURCES AND THE 
y-RAY log W-log S RELATION 

In Figure 1, we plot the number of blazars observed per 
square degree versus blazar flux integrated above 100 MeV for 
both EGRET (Reimer & Thompson 2001) and Fermi (Abdo 
et al. 201 Of). In the case of Fermi, the source spectra were 
extrapolated from a power-law fit above a fiducial energy of 
1 GeV (Abdo et al. 2010k). The offset between the resolved 
source count data obtained by the two detectors is a result of 
differences in their sensitivity calibrations (Thompson et al. 
1993; Abdo et al. 2010f). The fluxes of the Fermi sources are 
reported to be systematically brighter than those reported by 
EGRET. 

We note that the 50% completeness in the Fermi survey is 
estimated to be reached at an effective flux limit of ~2 x 
10 -8 cm -2 s -1 (Abdo et al. 2010k). The flux limit varies 
with galactic latitude and source hardness. Nevertheless, the 
effective flux limit indicated by both the turnover in Figure 1 
and the Monte Carlo simulation in Abdo et al. (2010k) is 
~2 x 10 -8 cm -2 s -1 . By assuming that the Fermi survey is 
100% complete at ~7 x 10 -8 cm -2 s -1 (Abdo et al. 2010k) 
and comparing the data for EGRET and Fermi shown in 
Figure 1, we estimate that the EGRET survey is 50% complete 
at ~8 x 10 -8 cm -2 s -1 . 

One may note that the flux level attained for 50% complete- 
ness for the Fermi extragalactic blazar survey is only about four 
times fainter than that for the EGRET survey even though the 
Fermi- LAT is sensitive to sources ~30 times fainter than that 
of EGRET. It is important to note that there are several factors 
that determine the efficiency of a y-ray telescope for detect- 
ing extragalactic sources, among which are (1) die flux of the 
source, (2) the spectral index of the source, (3) the intrinsic 
detector background from cosmic-ray-induced events, (4) the 
foreground from the Milky Way, and (5) the diffuse extragalactic 
background. The Fermi- LAT was designed to reach its optimal 


effective area for y-rays with energies near and above 1 GeV, 
whereas EGRET was designed to reach its optimal effective 
area for y-rays with energies near and above 100 MeV. As such, 
Fermi is more sensitive to sources with hard spectral indices, 
particularly for the faintest sources observable by Fermi. Thus, 
it is difficult to make a direct comparison between EGRET and 
Fermi. However, we note that the positions of the turnovers in 
the data presented in Figure 1 provide a good indication of the 
observational situation. 

3. DETERMINATION OF THE EGB FROM POINT 
SOURCES 

In general, the total contribution of a given population of 
sources to the EGB is found by convolving the source spectrum, 
F p h(Eo, z, L y ), measured at energy, E 0 , of one source of 
y-ray luminosity, L y , at a given redshift, z, with the comoving 
number density of sources at that redshift per luminosity 
interval, n(L y ,z) = d 2 N/dL y dV com . We then integrate over 
the comoving volume, V con „ and y-ray luminosity; 

fZma rZy.ma 

Ie(E 0 )= / / Fp h (Eo,z,L y )n<iL y ,z)- r ~dL Y dz, 

Jo Ji v dzdLl 

( 1 ) 

where Ie(Eq) is the EGB intensity given in units of photons 
cm -2 s -1 sr -1 GeV -1 , L Yim!a depends on the redshift, source 
spectrum, and the detector sensitivity, and we have differentiated 
the comoving volume with respect to redshift and solid angle, 
Q. However, as the different types of likely contributing y-ray 
sources have distinct spectral characteristics and redshift and 
luminosity distributions, we must consider each case separately. 

Note that in Equation (1) we have neglected the effect of 
y-ray absorption due to pair-production interactions with the 
extragalactic UV photons. The inclusion of y-ray absorption 
would result in a steepening in the collective point source 
spectrum at the high-energy end (Salamon & Stecker 1998; 
Venters et al. 2009), though a possible contribution from 
electromagnetic cascade photons might mitigate the steepening 
at higher energies (Venters 2010). However, recent Fermi 
constraints on the y-ray opacity imply that the UV background 
is likely to be fairly low (Abdo et al. 2010e), and as such, the 
absorption and the resulting cascades will only have a small 
effect on the EGB. 

3.1. The Contribution to the EGB from Unresolved Blazars 

In determining the blazar contribution to the EGB, we follow 
the procedure outlined in Venters et al. (2009). We approximate 
the blazar y-ray spectrum as a power law in energy defined 
by the photon spectral index, T (F ph « £ -r ). The spectral 
indices of the population of blazars form a distribution with 
some spread (the spectral index distribution, SID; Stecker & 
Salamon 1996; Venters & Pavlidou 2007); hence, the number 
density of blazars is defined as 

d 3 N 

n{L y ,z, T) = p y (L y , z)pl(T) = , (2) 

uLy d Vcom^ I 

where p y (L y ,z) — d 2 N /dL y dV (Xm is the blazar y-ray lumi- 
nosity function (GLF) giving the comoving number density of 
blazars per luminosity interval and Pl(V) = dN/dV is the nor- 
malized blazar SID accounting for spectral bias (see Section 4). 
L y is the y-ray luminosity at the fiducial energy, Ef (taken to 
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be 100 MeV), defined as Ej times the differential photon lumi- 
nosity, L p i = d 2 N Y /dtdE measured at £/. L y is related to the 
integral flux greater than Ef, £(>£/), by 

L y =47r£> 2 (T-l)(l + zf £/£(>£/), (3) 


where D is the distance measure for the Friedman-Robertson- 
Walker cosmology: 6 

L> = 4r f [Oa +« m (l +z') 3 ]~ 1/2 dz'. (4) 

Ho JO 

A given blazar of y-ray luminosity, L y , at a redshift, z, with a 
y- ray photon spectral index, T, has a measured photon flux of 


£ph(£o, z, Ly, r) — 


Ly 

4 rr£ 2 [<4(z)] 2 


(l + z) 2 - r 



, (5) 


where di{ z) = D( 1 + z) is the luminosity distance. Thus, 
the total contribution to the EGB at a given energy, £o, from 
unresolved blazars (the collective unresolved blazar intensity) is 
determined by integrating the contribution from each individual 
blazar fainter than the detector sensitivity, 


r oo rZaxx l* Ly muY 

/ 61 (^o)= / / / F^{Eo,Z,Ly,T) P y PL (T) 

J — OO *0 J l>yjn in ' 

( 6 ) 

where z m ax = 5.0 and Ly iimx is determined from Equation (3) 
taking £(>£/) — £mi n , with being the minimum flux 
capable of being resolved by Fermi. We have also integrated 
Equation (1) over the blazar spectral index. 

3.2. The Contribution to the EGB from Unresolved 
Star-forming Galaxies 

By applying the same procedure that we use to determine 
the background from unresolved blazars (see Section 3.1), we 
calculate the contribution from unresolved star-forming galaxies 
by determining the y-ray photon flux as.a function of energy 
for one galaxy and then convolving with and integrating over 
the appropriate cosmological distributions. It is expected that, 
as in the Milky Way, the y-ray emission for a star-forming 
galaxy comes mainly from the decay of n° mesons produced 
by cosmic-ray interactions with interstellar gas (Stecker 1970). 
The resulting y-ray production spectrum has been calculated 
by many authors (Stecker 1970, 1973, 1979; Cavallo & Gould 
1971; Stephens & Badhwar 1981; Dermer 1986; Mori 1997; 
Strong et al. 2004b, 2007, 2010; Kelner et al. 2006; Kamae 
et al. 2006; Mori 2009). For our calculation, we adopt the jt° 
emissivity given by Stecker (1979) renormalized upward by 
25% to be consistent with the local emissivity measured by 
Fermi (Abdo et al. 2009b). We note that there is also emission 
arising from electron bremsstrahlung, but the contribution is 
likely to be small, particularly above 100 MeV (Abdo et al. 
2009b). We also neglect the emission from Compton interactions 
that may contribute significantly to galaxy spectra above 10 
GeV, particularly for starburst galaxies (Stecker 1977; Hunter 
et al. 1997; Abdo et al. 2009b; Strong et al. 2010). We also 
note that the cosmic-ray spectrum should, in fact, vary from 


galaxy to galaxy since it depends on energy-dependent leakage 
and each galaxy has a different morphology and magnetic field 
configuration. This uncertainty can affect the predicted slope of 
the EGB spectrum at high energies, but it does not affect the 
absolute value of the predicted background at ~200MeV. 

The y-ray photon luminosity is related to the y-ray produc- 
tion spectrum per hydrogen atom, qn(E), by 

£ P h(£) = <9 h(£)) N h , (7) 

where (qu(E)) is found by averaging g H (£), the differential 
y-ray production spectrum per hydrogen atom, over the galaxy, 
and Nh is the number of hydrogen atoms in the galaxy in the 
form of both atomic (H i) and molecular (H 2 ) hydrogen. Thus, 
the y-ray photon flux for a galaxy at redshift z as observed at 
energy £0 is 

£ P h(£o. z) = («h[£o( 1 + z)]) Mt(z) • (8) 

The rate of production of y-rays from 7r° decay is proportional 
to the flux of cosmic rays, which we assume to be proportional 
to the supernova rate. The supernova rate is expected to be 
proportional to the rate of formation of higher mass stars, which 
is, in turn, proportional to the overall star formation rate (SFR) 
assuming a universal initial mass function (IMF). 

Assuming then that the rate of production of y-rays from jr° 
decay is proportional to the SFR for a galaxy, we can relate the 
(qn) of the galaxy to that of the Milky Way, <9 h W ): 

(«D v<z=or w 


where 4* is the SFR of the galaxy, and we take to be 

some fraction,/,, of the locally measured 7 q^- We determine 
f q by integrating the radial profile of the flux of cosmic rays 
weighted by r 2 . Using the radial profiles calculated by Stecker & 
Jones (1977), we obtain /, ~ 0.825. Thus, the galaxy spectrum 
becomes 

Fph(£ °’ z> = 4 ^ /9 ^ [£ ° (i + . go) 

The calculation of the amount of gas in a galaxy is subject 
to a considerable degree of uncertainty, especially at high 
redshifts (see Section 4.2). As such, rather than focusing on 
one particular model, we calculate the star-forming galaxy 
contribution for three different models arising from different 
sets of assumptions. In so doing, we seek to explore various 
possibilities and highlight the uncertainty. 

3.2.1. Galaxy Contribution Determined from the Schechter 
Function and An Evolving Gas Fraction 

One method for determining the number of hydrogen atoms 
in a galaxy is to assume that the mass of gas in the galaxy is 
some fraction of its stellar mass: 


V H (z) = 


/gas(z) M* 

1 “ /gas(z) mu 


(ID 


where / g2S (z) = M g aj/(Afgas + M.) oc (1 + z) 0 - 9 is the gas 
fraction given in Papovich et al. (201 1), and we have neglected 


4 We take //q = 70 km s 1 Mpc 1 , Q m = 0.3, = 0-7, and 0,4:1. 
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the possible contribution of helium to the gas mass of a 
galaxy. Substituting Equation (11) into Equation (10) and 
approximating 'P(z)/'F(0) ~ Psfr(z)/Psfr(0), we get 


F pb (£o. z ■ Af*) 


/,^ W [£o(l+z)] Psfr(z) / gM (z) 

4 nm H D 2 Psfr(0) 1 - / gas (z) 


( 12 ) 

where psfr(z) is the cosmic SFR (CSFR) given by 
log(/0SFR(z)) = -2.06 + 3.391og(l + z) for z < 1.3 and 
Psfr(z) ~ const, for 1.3 < z < 4.0 (Ly et al. 201 1). 

To get the total contribution to the EGB, we convolve with 
the comoving number density of star-forming galaxies per stellar 
mass interval as a function of redshift and integrate 


if (Eo) = 


d Z V, corn 
dzdCl 

/•"U 

x / dM' F ph (E 0 , z, M'mz, M') . 

J Kla 



(13) 


where ©(z, M') = d 1 N/dM'dV aom is the Schechter function 
for stellar mass with parameters as determined in Eisner et al. 
(2008), M' is given by M * = \0 M ' Mq, and we take M' m -„ — 8.0 
and = 12.0. 


3.2.2. Galaxy Contribution Determined from IR Luminosity Functions 


Alternatively, we can determine the y- ray spectrum of a 
galaxy by assuming that the y-ray luminosity of the galaxy 
is proportional to some power of its SFR (Fields et al. 2010; 
Makiya et al. 2011; Abdo et al. 2010g): L P h a 'F“. Since 
i, ph (£) = (q u (E)) Nh and (^h(^)> ex 'F (as demonstrated in 
Section 3.2.1), 

<14) 

where A and a are the best-fit parameters of the above power law 
determined from Fermi observations of star-forming galaxies in 
the Local Group and their SFRs (see Section 4.2.2). Assuming 
the Chabrier (2003) IMF, the SFR of a galaxy is related to its 
total infrared luminosity, Ljr, 

= 1.1X10 <15) 

(Hopkins et al. 2010). The y- ray photon flux for a galaxy at 
redshift z it given by 

Fph(£ °’ *’ Lir) = 4nD* 

Ku^) %n£o(1+z)] - 

(16) 


Then, the total contribution to the EGB is found by convolving 
the galaxy photon flux with an infrared luminosity function, 
©(z, Ljr) = d 2 N /d L^dVcotn and integrating over infrared 
luminosity and redshift: 


ifm = 



d 2 Vcom 

dzd£l 


x 


L 


^IR.nu? 

^IR.aiin 


dL m Fph(£o, z, Ljr)0(z, Lir), 


(17) 


where we take L^min — 10 10 L© and LiR.max — 10 15 L©. 


3.2.3. Galaxy Contribution Determined from the Cosmic Star 
Formation Rate and the Star Formation Efficiency 

Another alternative is to relate the cosmic density of hydrogen 
in star-forming galaxies to the cosmic SFR. Given that stars are 
formed in giant molecular clouds (GMCs), it is reasonable to 
assume 

Psfr ~ f (H 2 )Ph 2 , (18) 

where psfr is the cosmic SFR density (see Section 3.2.1), 
Ph 2 is the cosmic molecular hydrogen density in star-forming 
galaxies, and ?(H 2 ) is the star formation “efficiency” (SFE) of 
molecular hydrogen (Bigiel et al. 2008; Gnedin et al. 2009; 
Bauermeister et al. 2010). Leroy et al. (2008) measure the SFE 
to be ~(5.25 ± 2.5) x 10~ 10 yr _1 and to be roughly constant 
over a wide range of conditions. 8 We can relate the density of 
atomic hydrogen to the density of molecular hydrogen through 
the average mass ratio of atomic and molecular hydrogen 
(7 Z = {M m /M H2 )) in star-forming galaxies, p H , ~ TZpn 2 ■ The 
average mass ratio of atomic and molecular hydrogen can be 
found by integrating radial profiles of the gas surface densities 
of star-forming galaxies found in Leroy et al. (2008), resulting 
in TZ ~ 0.9. Note that in doing so, we only integrate the 
profiles out to the optical radius since recent surveys indicate 
that star formation is extremely inefficient beyond this radius 
(Bigiel et al. 2010). 9 Thus, with appropriate modifications to 
Equation (10), we find the y-ray flux from a particular redshift 


F p h(Eo, z) = 


fqd+K) 

4nm H |(H 2 )D 2 


s’TtEod+z)] 


Asfr(z) d Vcom 
Psfr(0) dz 


dz- 


(19) 

Differentiating with respect to solid angle Q, and integrating 
over redshift results in an equation for the total contribution to 
the EGB: 


E 4n (H 2 )psfr(0) Jo 

x + z)]p s 2 F R (z)^f^z. (20) 

4. OBSERVATIONAL INPUTS AND CONSIDERATIONS 
4.1. y-ray Blazars 

Of the nearly 1500 resolved point sources observed by 
Fermi in the first year, 573 are associated with blazars (First 
Fermi catalog (1FGL); Abdo et al. 2010f). Thus, blazars com- 
prise the largest class of astrophysical objects associated with 
y-ray sources. Naturally, unresolved blazars have long been 
suspected of providing, at least, a substantial contribution to the 
EGB, though the exact amount remains debatable and depends 
on various assumptions as to constructing GLFs and redshift dis- 
tributions (Padovani et al. 1993; Stecker et al. 1993; Salamon 
& Stecker 1994; Chiang et al. 1995; Stecker & Salamon 1996; 


8 More precisely, £(H 2 ) ~ e/rj, where € is the percentage of gas involved in 
forming stars and zs oc p -1 / 2 is the local free-fall timescale of the gas. 
However, using this relation requires knowledge of the local density of the gas 
and a better understanding of the formation of GMCs than presently exists (see 
Section 4.2). We also note that the measurements obtained by Leroy et al. 
(2008) were taken from a sample of low-redshift galaxies. The SFE could 
actually evolve with redshift (Bauermeister et al. 2010). 

9 We should note that even though star formation is extremely inefficient 
beyond the isophotal radius, there is still gas beyond this radius. However, ' 
based on the radial profiles determined in Stecker & Jones (1977), we do not 
expect cosmic rays to propagate much beyond the isophotal radius; hence, we 
expect gamma-ray production beyond this radius to be low. 
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Figure 2. Bright end of the source count distributions for all blazars (including 
BL Lac objects; black data points) and FSRQs (light green data points). Our 
model fit to the data is shown by the solid (purple) line. The dashed lines are 
the faint-end slopes determined in Abdo et al. (2010k) by including a modeled 
Monte Carlo Fermi-LAT efficiency. 

(A color version of this figure is available in the online journal.) 


Kazanas & Perlman 1997; Chiang & Mukherjee 1998; Sreeku- 
mar et al. 1998; Mukhetjee & Chiang 1999; Miicke & Pohl 
2000; Giommi et al. 2006; Narumoto & Totani 2006; Denner 
2007; Kneiske & Mannheim 2008; Pavlidou & Venters 2008; 
Inoue & Totani 2009; Venters et al. 2009; Abdo et al. 2010k; 
Venters 2010). A detailed discussion of all of the assumptions 
that go into these calculations is beyond the scope of this paper 
(though, for a detailed discussion of the Chiang & Mukheijee 
1998 calculation, see Stecker & Salamon 2001), and it is likely 
that many of these models will be updated in light of Fermi data. 
However, we note that as discussed in Abdo et al. (2010k), the 
source counts predicted by Dermer (2007) and Miicke & Pohl 
(2000) fall short of the Fermi observations of resolved sources 
above 5 x 10 -8 photons cm -2 s -1 . In calculating the blazar con- 
tribution to the EGB, we assume functional forms for the blazar 
GLF and SID and fit them to 1FGL data, accounting for errors 
in measurement of blazar spectral indices and the spectral bias 
inherent in a flux-limited catalog. 

4.1.1. Source Counts for Faint Unresolved Blazars: 

Theory Meets Observations 

Determining the GLF from observations relies on the ability 
to associate y- ray blazars with lower energy counterparts for 
which redshifts can be measured (for a discussion, see Venters 
et al. 2009; Venters 2010). However, making the necessary 
association can be complicated by the angular resolution of 
the Fermi-LAT, which is limited by electron scattering in the 
LAT detector and is much poorer than that of more traditional 
telescopes (see Figure 3). The resulting wide point-spread 
function results in significant source confusion at energies below 
~1 GeV even for fluxes well above the Fermi-LAT sensitivity. 
In principle, one could construct source counts from fluxes 
integrated above an energy for which source confusion is less 



Figure 3. Angular resolution as a function of energy for the Fermi-LAT (solid 
black line; Atwood et al. 2009) and EGRET (dashed red line; Thompson et al. 
1993). 

(A color version of this figure is available in the online journal.) 

of a hindrance, but doing so would limit the already suppressed 
blazar number statistics. Thus, rather than construct a luminosity 
function solely from y-ray blazars with redshifts, we employ 
the Stecker & Salamon (1996) 10 approach of determining the 
luminosity function from wavebands with larger samples and 
smaller positional error circles. 

As in Stecker & Salamon (1996), we take the functional 
form of the FSRQ luminosity function from radio observations 
(Dunlop & Peacock 1990), but corrected for the present cos- 
mological parameters. The y-ray luminosity of a blazar is then 
determined from its radio luminosity. The average correlation 
between the radio and y-ray luminosities of blazars is deter- 
mined by fitting the bright end of modeled source counts 11 to that 
of the observed y-ray source counts (x 2 duced ~ 0.4). In so do- 
ing, we find that the y-ray luminosity integrated from 100 MeV 
to 100 GeV is ~10 3 2 times vL v in radio, in agreement with 
the results obtained using recent Fermi observations 12 (Giroletti 
et al. 2010; Abdo et al. 2010j; Ghirlanda et al. 2011; Mahony 
et al. 2010). We should also note that we include sources out to 
z ~ 5, but as demonstrated in Venters et al. (2009), the emission 
is dominated by sources with z < 2 (consistent with expecta- 
tions from the redshift distribution predicted by the GLF). As 


10 See also Narumoto & Totani (2006). 

1 1 In doing so, we also include the blazar SID (see Section 4.1,3). 

12 One might be concerned that the effect of the new radio-y correlation 
would be to increase the blazar background with respect to the Stecker & 
Salamon (1996) model. However, we note that Stecker & Salamon (1996) 
distinguished between “quiescent” and “flaring” blazars and used separate 
radio-y correlations (and spectral properties) for each subpopulation. In this 
paper, we make no such distinction, so a comparison between our results and 
those of the Stecker & Salamon (1996) model is not straightforward. Using the 
Stecker & Salamon (1996) radio-y correlation for quiescent blazars from 
Stecker & Salamon (1996) would yield a smaller blazar background, but it 
would also underpredict the bright end of the y-ray source counts. To 
compensate, one would have to add a flaring component (as per Stecker & 
Salamon 1996) to fit the data, which would also contribute to the background. 
In effect, such a procedure is equivalent to our method of fitting the radio-y 
correlation to the data. 
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such, we do not expect this choice to significantly impact the 
results. 

The resulting modeled source count distribution for FSRQs 
(solid) is presented in Figure 2 along with the distributions of 
the brighter FSRQs resolved by Fermi (light data points) and 
all blazars (dark data points). We also show the forms of the 
unresolved source count distributions obtained from a Monte 
Carlo modeled Fermi-LAT sensitivity calculation performed by 
Abdo et al. (2010k) as dashed lines. All of the models separate 
from the data at fainter fluxes as the source counts fall off very 
rapidly due to the survey incompleteness and source confusion. 
Thus, the determination of the true faint-end shape directly from 
the source counts can be severely hindered. 

In an effort to mitigate the effect of survey completeness, 
Abdo et al (2010k) modeled the Fermi-LAT efficiency from 
a Monte Carlo simulation and then divided the differential 
source counts by this efficiency. Thus, the calculated source 
counts for faint sources strongly depend on such modeled 
efficiencies. Indeed, at fluxes near the Fermi-LAT sensitivity 
limit, the efficiency is extremely small and model dependent; 
hence, in flux bins with low number statistics, the source counts 
are multiplied by a very large and uncertain number. The result 
of such a procedure is that even though one source is not 
statistically different from two sources, whether one source 
is seen or two could result in different modeled counts for 
faint sources. It has also been argued that unassociated y- ray 
sources are likely to be dominated by known classes of y-ray 
sources, especially blazars (Mirabal et al. 2010) and will have 
a contribution to the EGB even though they will not have been 
included in the source count distributions. It is also important to 
note that in determining the blazar contribution to the EGB 
from measured source count distributions, source confusion 
could not be taken into account since its effect depends on 
the source density (see Section 4.1.2), which is exactly the 
unknown quantity that the observer seeks to determine. As 
such, it is likely that the blazar contribution to the EGB will be 
underestimated in analyses based solely on the measured source 
count distributions. Since in this paper we use a theoretically 
determined source count distribution, our model gives a source 
count density from which one can determine the effect of source 
confusion (see Section 4.1.2). Such differences between our 
analysis and that given in Abdo et al. (2010k) result in different 
calculations of the unresolved blazar contribution to the EGB: 
Abdo et al. (2010k) conclude that blazars can only account for 
less than 25% of the EGB, 13 while our analysis indicates that 
blazars could possibly account for the bulk of the EGB. 

The Fioc fluxes included in the source count distributions 
presented in Figure 1 are not actually measured Fioo fluxes 
since the 1FGL catalog does not include Fioo fluxes. In order 
to determine the Fioo fluxes, we make use of Equation (1) of 
Abdo et al. (2010k) to extrapolate Fioo fluxes from measured 
differential fluxes determined at the pivot energies. Abdo et al. 
(2010f) define the pivot energy as that energy for which the 
differential flux is minimal. We note that for 1FGL blazars, the 
average of the pivot energies is ~1 GeV. Thus, the source count 
distributions might be more representative of source counts for 
sources brighter than the Fermi sensitivity above 1 GeV rather 
than lOOMeV. In effect (and also as the result of spectral bias; 
see Section 4.1.3), the source counts could underestimate the 
number of blazars with Fioo fluxes above the Fermi sensitivity. 


mostly impacting the faint end of the source count distributions. 
As such, analyses on unresolved blazars based solely on source 
counts likely underestimate the blazar contribution to the EGB. 
This effect could provide an explanation for the fact that Fermi 
observes roughly as many BL Lac objects as FSRQs even though 
BL Lac objects are intrinsically fainter than FSRQs. Since the 
y- ray spectra of BL Lac objects are harder than those of FSRQs 
they are more easily observed by Fermi than FSRQs. 

4.1.2. The Fermi-LAT and EGRET Angular Resolutions 
and Source Confusion 

As mentioned previously, the large angular resolutions of 
pair-production y-ray detectors such as EGRET and the Fermi- 
LAT 14 result in significant source confusion, particularly for 
faint sources and for energies below ~1 GeV. Since the angular 
resolution for the Fermi- LAT is similar to that of EGRET at 
'-“100-200 MeV, the capability of Fermi to resolve faint sources 
is similar to that of EGRET at these energies. Hence, if the 
EGB does indeed consist of unresolved sources, then the EGB 
measurements of the detectors should be similar at ~100- 
200 MeV, whereas at ~1 GeV, the improved angular resolution 
of Fermi with respect to that of EGRET should result in a lower 
measurement of the EGB by Fermi than that of EGRET owing to 
the enhanced ability of Fermi to resolve point sources (Stecker 
& Salamon 1999). 

The probability, V, of finding a nearest neighboring source 
with S ^ Sum within the minimum angular separation, 0 m i n , for 
a source density, N, is given by 

flmin (£)) = 1 - exp (-jtV^(F)). (21) 

For our source confusion criterion, we take the acceptable prob- 
ability limit, Vaun ~ 0.1, and 0min(F) ~ 0 61 %(E) approximately 
given by 


<22> 

(Atwood et al. 2009; see Figure 3). Then, the source density 
criterion (SDC) is found by inverting Equation (21): 


IVsdc - — 


ln(l - Pnun) 

X8LW ' 


(23) 


The limiting source flux, Ssdc> is then determined from the 
modeled source counts. 

If one were to think of the Fermi-LAT as an ordinary tele- 
scope with 0mi n = 067 % where 067 % is the half-angle for 
a beam that contains 67% of the photons, the SDC would 
correspond to M/10 sources per beam. For sources with 
Fioo ;$ 1 x 10 -7 photons cm" 2 s -1 , the probability for find- 
ing another y- ray source of similar or greater flux within the 
error circle is quite high; hence, many sources that, in prin- 
ciple, should be resolvable are, in fact, unresolved. At fluxes 
close to the sensitivity limit, this probability is so high that 
faint sources are indistinguishable and will contribute to the 
measured EGB of the detector in question (for the Fermi- 
LAT resolution at lOOMeV, the source criterion corresponds 
to Nsdc ~ L3 x 10 -3 sources deg -2 ; for reference, at F ~ 
2 x 10 -9 photons cm -2 s -1 , our model predicts a much larger 
source density of ~0.3 sources deg -2 ). Thus, the effect of source 


13 In Section 8 of Abdo et al. (2010k) a value of 40% of the EGB is obtained if 14 See http://www-glast.slac.stanford.edu/soflware/IS/glast_lat_ 

one extrapolates to zero flux. performance.htm. 
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confusion compounded with the separate effect of detector sen- 
sitivity is to flatten the faint end of an observationaUy derived 
source count distribution. 

We should emphasize that this definition of source confusion 
is not the same as that employed by Abdo et al. (2010k). Source 
confusion as discussed in this section refers to the probability 
that for a given source density, a source has a nearest neighbor 
within the angular resolution of the detector with a flux greater 
than or equal to the flux limit. In Abdo et al. (2010k), the 
term is applied to a detected source associated with a given 
real source for which the measured flux of the detected source 
is greater than the flux of the real source (plus three standard 
deviations) by a given amount (the analysis was performed on 
simulated data to determine the impact on the actual data). In 
effect, the former definition identifies the limit at which the 
source density is sufficiently large so that individual sources 
cannot be resolved, and only fluctuations are observed. The 
latter definition applies to the probability that a given detected 
source is in fact the superposition of several sources. However, 
the Abdo et al. (2010k) criterion is based on the assumption 
that a source can be resolved. As such, the source would have 
to be significantly brighter than the background, including the 
sources within the error circle of the detector. In the case when 
many sources have similar fluxes, none of the sources would 
be resolved. In the case of one bright source and several fainter 
sources, the flux of the detected source would be dominated 
by the flux of the brightest source, leading to an overestimate 
of the flux of the dominant source and an undercount of the 
fainter source in the Abdo et al. (2010k) analysis. Thus, under 
this criterion, it is not surprising that they conclude that there 
is very little source confusion. A common treatment of source 
confusion in measured source counts is a fluctuation analysis, 
but as yet, such an analysis has not been performed on Fermi 
source counts. 


4.1.3. Blazar Spectra 

The spectral indices of the population of blazars form a 
distribution with a given finite spread (Stecker & Salamon 
1996; Venters & Pavlidou 2007), resulting in a curvature in the 
spectrum of the collective intensity of unresolved blazars due 
to the increasing relative importance at high energies of blazars 
with harder spectral indices (Stecker & Salamon 1996; Pavlidou 
& Venters 2008). Thus, the determination of the blazar SID is 
crucial in determining the correct spectrum of the unresolved 
blazar contribution. 

In determining the blazar SID from survey data, one must 
carefully account for uncertainties in measurement of the 
spectral indices. In so doing, we follow the likelihood method 
of Venters & Pavlidou (2007), fitting Fermi 1FGL FSRQs 
to a Gaussian SID. We determined the maximum-likelihood 
Gaussian SID parameters (mean. To, and spread, oo) to be 
To = 2.45 and cro = 0.15. 

We must also account for the effect of spectral bias of the 
1FGL catalog. Even in a flux-limited catalog, low-luminosity, 
high-redshift blazars that would be most likely to appear are 
those with spectral indices that are harder than most of the 
population. The 1FGL catalog is not actually flux limited (see 
Section 2); rather, it consists of sources above a given threshold 
in test statistic, 15 which depends on source spectra and the 
background. However, as demonstrated in Venters & Pavlidou 


15 For 1FGL, this threshold is 25, corresponding to a statistical significance of 
^4 <T. 


(2010), applying the likelihood analysis to the sample of FSRQs 
with fluxes >7 x 10~ 8 photons cm -2 s _l and galactic latitudes 
>10° (as per Abdo et al. 2010k) does not appreciably change the 
SID parameters (To = 2.45, cro = 0.16). 16 Thus, even though 
the 1FGL catalog is not flux limited, we can assume that the 
sample of 1FGL FSRQs is approximately flux limited. We can 
then follow the method of Venters et al. (2009) in correcting for 
the sample bias inherent in a flux-limited catalog. In doing so, 
we apply a correction factor, Mia), to the SID (for derivation, 
see Venters et al. 2009): 


Pi. (a) = 


P(«) 
M{a) ’ 


(24) 


where /3(a) is the SID corrected for measurement uncertainty 
in the spectral indices, and 


M{a) oc f dF y ~ f dzp y (a,z, F y )^~^{z). (25) 
r y Jz=a dz 

4.1.4. Summary of Differences with the Stecker & Salamon 
( 1996) Blazar Model 

To summarize, our approach in calculating the blazar con- 
tribution to the EGB is similar to that of Stecker & Salamon 
(1996) with some notable differences. 

1. The model has been updated to make use of the current 
cosmological parameters. The change in cosmology has 
little impact on the results. 

2. The model considers only FSRQs, and flaring blazars 
are not considered separately from quiescent blazars. (For 
possible impact, see item 3.) 

3. The y-ray-radio relation has been updated to be consistent 
with multi-wavelength observations of FSRQs conducted 
by Fermi and radio telescopes (L y ~ 10 3 - 2 x vL v as op- 
posed to 10 2 6 ; see Section 4. 1 . 1). As noted in Section 4.1.1, 
the smaller y-ray-radio relation would result in a model 
that cannot by itself fit the data and would require a flar- 
ing component. In effect, this is, on average, equivalent to 
our choice of a single population of blazars with a given 
y-ray-radio relation. We do not expect this choice to have 
a major impact on our results. However, as we will discuss 
in Section 4.1.5, the blazar duty cycle is a remaining un- 
certainty, and it could impact predictions for the number of 
blazars that are observable by Fermi. 

4. The model has been updated to account for the effect of 
source confusion in the blazar contribution to the EGB (see 
Sections 4.1.2 and 5). As will be shown in Section 5, this 
has the effect of increasing the blazar background at lower 
energies. 

5. The SED of FSRQs has been updated following the analysis 
of Venters & Pavlidou (2007; To = 2.45, a 0 = 0.15) and 
correcting for spectral bias as in Venters et al. (2009; see 
Section 4.1.3). Thus, the collective spectrum of blazars is 
not as hard as that presented in Stecker & Salamon (1996) 
and does not exhibit as much curvature. 

4.1.5. Remaining Questions 

There remain a few open questions, the answers to which will 
impact the determination of the blazar contribution to the EGB. 
The blazar duty cycle, which dictates the amount of time a blazar 


16 See also Abdo et al. (2010k). 
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spends in the quiescent state versus the flaring state, remains 
uncertain, as do questions of the amount the flux increases during 
flaring and whether the spectral index changes during flaring. 
Analyses of EGRET blazar spectral indices found no evidence of 
systematic changes in spectral index with flaring (Nandikotkur 
et al. 2007 ; Venters & Pavlidou 2007), and Fermi observations of 
individual blazars have thus far revealed no systematic changes 
in spectral index with time or flux (Abdo et al. 2009c, 2010d; 
Ackermann et al. 2010b). As such, we are justified in assuming 
that the blazar spectral index remains constant, on average, with 
time. However, we acknowledge that the uncertainty of blazar 
variability parameters could have an impact on the counts of 
faint blazars and the y -ray-radio correlation. This uncertainty 
will decrease as more data from Fermi become available. 

Another open question is that of the nature of blazar spectra 
over the entire Fermi energy range. We treat blazar spectra as 
unbroken power laws over this range and in many observed 
blazars, this does appear to be a reasonable approximation. 
However, in at least a few cases, Fermi has found evidence 
that blazar spectra can break (Abdo et al. 2009a, 2010i, 20101). 
Whether such observations are representative of the entire blazar 
population is presently unclear, as is the nature of the breaks. 
In any case, .spectral breaks are likely to impact the collective 
unresolved blazar spectrum mostly at the high end of the Fermi 
energy range. It is also possible that the spectra of harder 
blazars 17 will compensate for that of softer blazars at higher 
energies (Venters & Pavlidou 2010). 

4.2. Star-forming Galaxies 

As discussed in Section 3.2, the Milky Way is a source of 
substantial y-ray emission arising primarily from the decay 
of 7r° meson produced in inelastic collisions of cosmic rays 
with interstellar gas. Thus, it is expected that other star-forming 
galaxies emit y-rays through the same interactions and that 
unresolved star-forming galaxies could provide a substantial 
contribution to the EGB . Thus far, the Fermi- LAT Collaboration 
has reported detections of two nearby irregular galaxies (the 
Small Magellanic Cloud (SMC) and the Large Magellanic Cloud 
(LMC); Abdo et al. 2010c, 201 Oh), two starburst galaxies (M82 
and NGC253; Abdo et al. 2010b), and M31, a galaxy similar to 
our own (Abdo et al. 2010g). As such, whatever the contribution 
to the EGB from unresolved star-forming galaxies, it will not 
have changed substantially in the Fermi data with respect to 
EGRET data as Fermi has resolved only a handful of star- 
forming galaxies. 

The determination of the star-forming galaxy contribution 
relies on knowledge of the SFR of galaxies and their gas con- 
tent, both of which are subject to substantial observational and 
theoretical uncertainties. At relatively low redshifts (z < 1.5), 
nebular and forbidden emission lines (e.g., Ha, On, and Om) 
can be used to trace star formation in galaxies. 18 At higher 
redshifts (z ~ 1-5), the redshifted UV continuum is used. How- 
ever, both observational techniques are subject to uncertainties 
in dust extinction and the stellar IMF. Alternatively, in noting 
that UV radiation from young stars is absorbed by interstellar 
dust and reradiated in the infrared, measurements of the far- 
infrared (FIR) continuum can be used to trace star formation. 
Nevertheless, early-type galaxies can exhibit substantial FIR 


17 We do not include a possible contribution from BL Lac objects for lack of a 
comparable radio data set. 

18 For a review of observational techniques of measuring star formation, see 
Kennicutt (1998). 


emission possibly due to dust heating by older stars or AGNs. 
Furthermore, infrared measurements are hindered by emission 
from our own Galaxy. Given these factors, the large degree of 
scatter present in the measurements (about a factor of a few; Le 
Borgne et al. 2009; Ly et al. 201 1) of the CSFR density is not 
surprising. 

The gas content of galaxies is even more uncertain, particu- 
larly at high redshift. The amount of H 2 in a galaxy is determined 
from measurements of CO emission, while Hi is determined 
from measurements of the 21 cm line. However, the CO-to-H 2 
conversion varies depending on the metallicity and radiation 
field of a given region and the opacity of the molecular clouds 
containing CO. The 21 cm surveys, on the other hand, extend 
only out to z ~ 0.05. At higher redshifts, measurements of the 
Hi density of the universe rely on damped Lya absorbers ob- 
served in the Lya forest of quasar spectra, but the nature of these 
systems is still the subject of much debate (see, e.g., Kulkami 
et al. 2010; Peroux et al. 2011). Furthermore, the connection 
between the total gas and SFR is complex (Putman et al. 2009). 
While it is fairly well established that stars form in GMCs and 
hence that star formation traces H 2 , the amount of Hi varies 
between galaxies and does not appear to be correlated with star 
formation (Bigiel et al. 2008; Leroy et al. 2008). From both the 
observational and the theoretical points of view, the transition 
from H 1 to H 2 and the formation of GMCs remain uncertain 
(Leroy et al. 2008). 

Given the uncertainty surrounding key elements of the deter- 
mination of the star-forming galaxy contribution to the EGB, 
our approach does not focus on a particular model. Instead, we 
employ several families of models that rely on separate sets 
of assumptions each with advantages and caveats. Using this 
approach, we seek to explore various possibilities for the star- 
forming galaxy contribution and highlight the uncertainty in 
such a calculation. The strategies we employ are summarized as 
follows. 

1. Relate the galaxy gas mass to its stellar mass assuming a 
gas fraction that evolves with redshift. 

2. Relate the galaxy y-ray luminosity to its SFR, which, in 
turn, is related to an observable for which there is a redshift 
distribution (e.g., IR luminosity). 

3. Relate the cosmic density of gas in star-forming galaxies to 
the SFR density. 

4.2. 1. The Schechter Function Model 

In this approach, we relate the galaxy gas mass to its stellar 
mass assuming an evolving gas fraction (for details of the 
model, see Section 3.2. 1). We employ the Schechter parameters 
of the stellar mass functions as determined by Eisner et al. 
(2008) and the evolving gas fraction as determined by Papovich 
et al. (201 1). Since extensive spectroscopic surveys are, as yet, 
unavailable, Eisner et al. (2008) make use of combined data from 
the multi-band photometry of the GOODS-MUSIC catalog and 
the Spitzer Space Telescope. In so doing, they infer the stellar 
masses of galaxies from photometric data by fitting the mass- 
to-light ratios of galaxies to stellar population templates. Such 
a procedure is subject to a considerable degree of uncertainty, 
particularly arising from that of dust extinction 19 and the usage 
of photometric redshifts. Eisner et al. (2008) estimate the mean 
uncertainty in their stellar mass estimates to be about a factor 


19 In order to mitigate the effect of the uncertainty in dust extinction, Eisner 
et al. (2008) make use of Xs-band M/L ratios since dust absorption is small for 
longer wavelengths. 
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Table 1 

Observables for Local Group Galaxies 


lD a 

L y (photons s 3 ) 

V(M a yr-‘) 

SMC 

(1.7 ±0.3) x IO 40 

0.01 ± (4.8 x IO” 4 ) 

LMC 

(7.8 ±0.6) x 10 40 

0.1 ± (4.6 x 10 -3 ) 

M31 

(6.6 ± 1.4) x 10 41 

0.4 ±0.18 

Milky Way 

(1.2 ±0.2) x 10 42 

0.93 ±0.34 

NGC 253 

(1.1 ±0.7) x 10 43 

2.7 ± 0.5 

M82 

(2.5 ±0.9) x 10 43 

5.4 ±(5.9 x IO" 4 ) 


Note. * SMC: L y from Abdo et al. (2010c). SFR from Lawton et al. (2010) 
corrected for IMF. LMC: L y from Abdo et al. (2010h). SFR from Lawton 
et al. (2010) corrected for IMF. M31: L y from Abdo et al. (2010g). SFR from 
Tabatabaei & Berkhuijseii (2010) and Williams (2003) corrected for IMF. MW: 
L y calculated for (qu) ~ 1.4 x IO -25 s _l H 1 and Mb ~ 7 x 10 9 M 0 (Boissier 
& Prantzos 1999). SFR taken from Robitaille & Whitney (2010) corrected for 
IMF. NGC 253: L y from Abdo et al. (2010b). SFR from 1R measurements 
given by Sanders et al. (2003). M82: L y from Abdo et al. (2010b). SFR from IR 
measurements given by Sanders et al. (2003). Note that the errors in the SFRs 
reflect only statistical uncertainties and assume that Equation (15) is exactly 
correct. 

of two, though they did not estimate the possible uncertainty 
resulting from the usage of photometric redshifts. 

Papovich et al. (2011) study the relationship between the 
SFR and stellar mass of high-redshift galaxies selected at a 
constant comoving number density and derive a gas fraction 
that evolves 20 as (1 + z) 0 9 . They estimate that the uncertainty 
on the gas mass is ~0. 1 1 dex. 

The advantage of this kind of model is that it is based 
on observations that will be continuously refined with time. 
However, we note that the relationship between the stellar mass 
of a galaxy and its total gas content is unclear, though perhaps 
with better observations and a better theoretical understanding, 
the relationship will be better determined. 

4.2.2. The IR Luminosity Function Models 

In this approach, we assume that the x-ray luminosity of a 
star-forming galaxy can be related to its SFR as a power law. In 
order to determine this relationship, we fit y- ray luminosities of 
Local Group galaxies calculated from Fermi measurements 21 to 
their SFRs either taken from the literature or calculated from IR 
measurements and converted to SFR via Equation (15). In all 
cases, the SFR is calculated assuming (or corrected to) the IMF 
of Chabrier (2003). We find the best-fit power law to be given by 
L y (10 41 photons s -1 ) ~ 24.0 x T 1,2 . This fit is consistent with 
that obtained by the Fermi Collaboration, L y <x 'p L4±0 - 3 (Abdo 
et al. 2010g). Observables for Local Group galaxies used in this 
analysis are given in Table 1 and plotted with the fit in Figure 4. 
In a manner s imil ar to that of blazars, we can relate the y- ray 
luminosity of a star-forming galaxy to its total IR luminosity, 
convolve with an IR luminosity function, and integrate with 
respect to IR luminosity and redshift. 

While the use of IR luminosity functions taken from obser- 
vations is possible, it is difficult to deconvolve the contribu- 
tions from obscured AGNs and mergers. As such, we use the 
semi-empirical IR luminosity functions determined from the 
halo-occupation-based methodology of Hopkins et al. (2010) 


20 Note that we extrapolate the functional form of their gas fraction to low 
redshifts. 

21 For the Milky Way, we calculate the y-ray luminosity from qyi averaged 
over the whole galaxy (as discussed in Section 3.2) and the Mb taken from 
Boissier & Prantzos (1999). 



Figure 4. y-ray luminosities of Local Group galaxies plotted vs. their star 
formation rates (data points; see Table 1). Also plotted is the power-law fit 
relating the y-ray luminosity to the SFR (solid red line): L y oc ¥ . 

(A color version of this figure is available in the online journal.) 

for both star-forming and starburst galaxies (in which enhanced 
star formation due to major mergers is taken into account). 22 
While the Hopkins et al. (2010) IR luminosity functions match 
available observations fairly well, we note that there is consid- 
erable debate over the roles of AGNs and mergers in driving 
star formation, the evolution of galaxies, and the determination 
of the IR luminosities of massive systems. Furthermore, the 
y-ray-SFR correlation is based on rather uncertain estimates 
of the SFRs and y-ray luminosities of one normal galaxy (our 
own), two irregular galaxies, and two starburst galaxies, all of 
which could, in principle, exhibit different star formation proper- 
ties. As more data become available from Fermi, the y-ray-SFR 
correlation will be further tested. If it proves robust, then studies 
of the normal galaxy contribution to the EGB could have impli- 
cations for large-scale-structure formation and the evolution of 
galaxies with cosmic time. 

4.2.3. The Strong Coupling y -ray-Star Formation Rate Model 

In this approach, we seek to relate the cosmic density 
of hydrogen in star-forming galaxies to the cosmic SFR. 
Observations of nearby galaxies have indicated that localized 
star formation traces the density of H 2 but there is no direct 
correlation with Hi (Bigiel et al. 2008; Leroy et al. 2008). 
While, in principle, there should be some relationship between 
the amount of Hi, the amount of H 2 , and the SFR in a galaxy, 
other factors such as density fluctuations and turbulence make. 
such a relationship complex. Therefore, while acknowledging 
that the relationship between Hi and SFR is unclear, we simply 
assume that the amount of H 1 in star-forming galaxies is, on 
average, comparable to that of H 2 within the optical radius 


22 Note that we do not consider any contribution from the so-called 
calorimetry effect for starburst galaxies as such an effect is likely to be small 
(Stecker 2007). We have assumed the same form of the jr°-decay spectrum for 
starburst galaxies as for star-forming galaxies. 
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Figure 5. Collective spectrum of unresolved FSRQs. Solid green line: the 
spectrum accounting for source confusion. Dashed green line: the spectrum 
without accounting for source confusion. Black circles: the Fermi measurement 
of the spectrum of the EGB as determined in Abdo et al. (2010j). Blue 
squares: the EGRET measurement of the spectrum of the EGB as determined 
by Sreekumar et al. (1998) and confirmed by the analysis of Stacker et al. 
(2008) and S. D. Hunter (2011, private communication). Red triangles: the 
EGRET measurement of the spectrum of the EGB as determined by Strong 
et al. (2004a). 

(A color version of this figure is available in the online journal.) 


of the galaxy (see Section 3.2.3). Hence, in this model, we 
take the y-ray luminosity to be roughly proportional to the 
square of the SFR. We stress that the assumption that the SFR is 
proportional to the available gas density only applies to galaxies 
that are actively forming stars. This assumption does not apply 
to galaxies at very high redshifts which may contain substantial 
amounts of gas but have yet to begin forming stars. However, 
we only include star-forming galaxies out to z ~ 4, so the 
higher redshift galaxies will not impact on our results. Given 
that the best-fit power-law index determined for Local Group 
galaxies in Section 4.2.2 is ~1 .2 and given the proximity of the 
resulting unresolved spectrum to the Fermi measurements of 
the EGB (see Section 5), we consider this model to be reflective 
of an upper limit to the star-forming galaxy contribution to 
the EGB. 

4.3. The Fermi Spectrum and Unresolved Sources Versus 
Truly Diffuse Mechanisms 

The Fermi observations have placed significant constraints 
on extragalactic dark matter annihilation (Cirelli et al. 2010; 
Abdo et al 2010a; Ackermann et al. 2010a). Currently, there 
is no evidence of quark-annihilation features and spectral lines 
seen in the EGB spectrum, features that would be a clear an- 
nihilation signal (see, e.g., Stecker & Tylka 1989a; Rudaz & 
Stecker 1991). The observed spectrum does not match that ex- 
pected front dark matter annihilation, placing constraints on 
any dark matter annihilation contribution to the EGB (Abdo 
et al. 2010a). Therefore, it is probable that dark matter annihi- 
lation y-ravs, if present, provide only a minor contribution to 
the EGB. 


The same argument about matching spectra can be made 
regarding the contribution from electromagnetic cascades pro- 
duced by very high and ultrahigh energy cosmic-ray interactions 
as the resulting spectrum would be significantly harder than the 
observed spectrum (Kalashev et al. 2009; Berezinsky et al. 201 1 ; 
Ahlers et al. 2010; Venters 2010). 

5. RESULTS 

The calculated spectrum of the unresolved FSRQ contri- 
bution to the EGB (see Sections 3.1 and 4.1) is plotted in 
Figure 5. For comparison, we include the Fermi analysis of the 
EGB (Abdo et al. 2010j), two analyses 23 of the EGRET EGB 
(Sreekumar et al. 1998; Strong et al. 2004a), and the calcula- 
tion of the collective spectrum of unresolved FSRQs ignoring 
the effect of source confusion. Our results clearly show that 
the effect of source confusion is to reduce the number of re- 
solved sources, increasing the collective intensity of unresolved 
blazars, particularly below ~1 GeV energy. Thus, accounting 
for source confusion modifies the predicted spectrum such that 
the EGRET and Fermi measurements of the EGB below ~1 
GeV are both compatible with unresolved FSRQs. In contrast, 
the better angular resolution of the Fermi-LAI above ~ 1 GeV 
allows it to resolve more blazars resulting in a limiting flux that 
is dominated by the Fermi-LAI sensitivity rather than source 
confusion. Thus, the collective spectrum of FSRQs breaks at 
~3 GeV. 24 At energies above ~1 GeV, the predicted collective 
spectrum of FSRQs falls below the data points, though they are 
likely consistent with the data within the uncertainties in the 
galactic foreground emission model. Note also that the collec- 
tive FSRQ spectrum exhibits much less curvature than seen in 
Stecker & Salamon (1996). This is because the spread in the SID 
in our current model is much smaller than that of the Stecker & 
Salamon (1996) model. 

In Figure 6, we plot the spectra of the unresolved star- 
forming galaxy contributions to the EGB calculated for the 
models discussed in Sections 3.2 and 4.2. For comparison, 
we include the spectrum of the unresolved starburst galaxy 
contribution alone that we determined from the best-fit IR 
luminosity function of Hopkins et al. (2010). For the spectrum 
of starburst galaxies, we have assumed the same form of the 
7T°-decay spectrum as for star-forming galaxies. The range 
in the calculations of the overall contribution to the EGB 
from unresolved star-forming galaxies spans about an order 
of magnitude, indicating the degree of uncertainty in such a 
calculation. 25 

We note that even though our most extreme model could 
possibly explain the lowest energy Fermi data points (and 
possibly, within systematics, a couple other), it cannot explain 
the EGRET data points below 300 MeV. The Strong et al. 
(2004a) EGRET data points (minus the two highest energy data 
points) with the Fermi data points resemble a . featureless power 
law, while the spectra of unresolved star-forming galaxies do 
not. Notably, the data points show no indication of a 7r°-decay 


23 The two sets of EGRET data points result from two different estimations of 
the galactic foreground emission. 

24 The actual break should be more gradual since in our calculations we used 
the approximate broken angular resolution curve shown in Figure 3. 

25 However, we note that each individual model is subject to its own 
uncertainty. As such, the degree of uncertainty is likely even more than an 
order of magnitude. Within the range of our various predictions of the EGB 
from star-forming galaxies, we agree with the results of the model of Fields 
et al. (2010) and Makiya et al. (201 1). 
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Figure 6. Collective spectrum of unresolved star-forming galaxies. Dashed 
indigo line: the spectrum determined from the strong coupling model (see 
Sections 3.2.3 and 4.2.3). Solid blue line: the spectrum determined from the 
IR luminosity function model (see Sections 3.2.2 and 4.2.2). Dot-dashed yellow 
line: the spectum determined from the IR luminosity function model assuming 
no gas evolution. Dashed red line: the spectrum determined from the Schechter 
function model (see Sections 3.2.1 and 4.2.1). Double dot-dashed line: the 
spectrum of tie starburst contribution alone determined from the IR luminosity 
function model. 

(A color version of this figure is available in the online journal.) 

“bump" at the energies at which the contribution of the star- 
forming galaxies should peak. 


6. DISCUSSION AND CONCLUSIONS 

We have calculated the spectral shape of the contribution 

of unresolved FSRQs to the EGB assuming that the y-ray 
luminosity of an FSRQ is, on average, proportional to its 
radio luminosity (Giroletti et al. 2010; Abdo et al. 2010j; 

Ghirlanda et al. 2011; Mahony et al. 2010), and also accounting 
for the effects of source confusion. We have demonstrated 

that the combination of the source density predicted by the 
Dunlop & Peacock (1990) FSRQ radio luminosity function 

and the strong energy dependence of the Fermi- LAT angular 
resolution increases the contribution of unresolved FSRQs to the 

EGB at energies below 1 GeV. The resulting overall spectrum 

predicted by the fit to the Fermi source count distribution 

reproduces well the spectrum of the EGRET and Fermi EGB 

measurements below 1 GeV, but falls below the data points 

above 1 GeV. We have also calculated the spectral shape of 

the contribution of unresolved star-forming galaxies to the 

EGB for several relations for the y-ray luminosity of a star- 

forming galaxy. We find that, depending on the model, the 

overall amount of the contribution of star-forming galaxies to 

the EGB may be more or less significant, though regardless 

of the model considered, the spectrum of unresolved star- 

forming galaxies is unable to explain the combined spectrum 

of the low-energy EGRET EGB measurements and the Fermi 

EGB measurements. Similar calculations for starburst galaxies 

alone indicate that they account for at most about 1% of the 

EGB, in agreement with the conclusion reached by Stecker 

(2007). 


The similarity of the collective spectrum of unresolved 
FSRQs to the combined spectrum of the EGRET and Fermi 
EGB measurements, as demonstrated by our results, is striking. 
In fact, we note that as predicted in Stecker & Salamon (1999), 
the inclusion of the effect of source confusion in the calculation 
could provide an explanation for the similarity between the 
EGRET and Fermi EGB measurements at energies of hundreds 
of MeV. The density of FSRQs predicted by the model is 
sufficiently large such that at these energies, Fermi would not 
be able to resolve many more FSRQs than EGRET did, and 
the FSRQ contribution to the EGB would remain the same for 
Fermi as for EGRET. Thus, if unresolved FSRQs do comprise 
the bulk of the EGB emission, then one would expect such 
a similarity between the EGRET and Fermi measurements at 
these energies. At energies above 1 GeV, the Ferroi-LAT angular 
resolution improves substantially with respect to thatofEGRET, 
and as such, Fermi would be able to resolve more blazars at 
higher energies than EGRET could, resulting in a decrease in 
the Fermi EGB with respect to the EGRET EGB, an effect 
which is possibly indicated by comparing the EGRET and Fermi 
results. 2 ® 

In contrast, no such high-energy separation between the 
EGRET EGB and the Fermi EGB is predicted for star-forming 
galaxies as all but the closest are too faint to be resolvable by 
Fermi. Furthermore, we note that the EGRET EGB measure- 
ments provide no indication of a turnover in the spectrum as 
would be expected if unresolved star-forming galaxies com- 
prise the bulk of the EGB emission. Rather, the spectrum of 
unresolved star-forming galaxies is inconsistent with the com- 
bined spectrum of the EGRET and Fermi EGB measurements. 6 * * * * * * * * * * * * * * * * * * * * 27 
We also note that the lack of a turnover in the EGRET data is 
not simply the result of systematics (S. D. Hunter 2011, private 
communication), since the uncertainties in all of the galactic 
foreground models used to determine the EGB from the EGRET 
and Fermi data are quite small at these energies. Finally, we note 
that at energies above ~1 GeV, the spectrum of unresolved star- 
forming galaxies is steeper than the spectra of the EGB data. As 
such, we conclude that however significant the contribution of 
star-forming galaxies to the EGB may be, it is not sufficient to 
explain the EGB. 28 

Within the range of our various predictions of the EGB from 
star-forming galaxies, we agree with the results of the models 
of both Fields et al. (2010; which suggests that star-forming 
galaxies may comprise the bulk of the EGB) and Makiya et al. 
(2011; which suggests that star-forming galaxies can account 


26 A caveat is that the uncertainties in the subtraction of the galactic 
foreground emission at the higher energies are considerable owing to the 
uncertainty in the distributions of both gas and cosmic rays in the Galaxy. 
Furthermore, the instrumental backgrounds of EGRET and the Fermi - LAT are 
different, so it is difficult to make a direct comparison between the two. We 
should also note that in our calculations, we have neglected the population of 
BL Lac objects, which, due to their hard spectra, are likely to have more of a 
contribution at energies above ~10 GeV. Notably, Fermi has resolved as many 
BL Lac objects as FSRQs. 

27 As previously noted, the Fermi- LAT was designed to reach its optimal 
effective area for x-rays with energies near and above ~1 GeV, whereas 
EGRET was designed to reach its optimal effective area for y-rays with 
energies near and above ~ 100 MeV. Also, the Fermi-LAT detector has a 
significantly higher instrumental background at 100 MeV than EGRET did 
(S. D. Hunter 2011, private communication). Thus, the EGB was not reported 
by Fermi for energies below 200 MeV (Abdo et al. 2010j). 

28 The effect of Compton interactions mentioned in Section 3.2 does not alter 
this conclusion as it only modifies the spectrum above 10 GeV for normal 
galaxies (Strong et ai. 2010) and the starburst galaxy contribution to the EGB 
is negligible. (See Figure 6.) 
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for less than 10% of the EGB). This underscores the range of 
uncertainty in the calculation for star-forming galaxies. 29 

The featureless spectrum of the EGB deduced by Fermi 
is intriguing when one considers the possibility of features 
that could arise from phenomena such as breaks in blazar 
spectra, absorption of high-energy y-rays from unresolved 
blazars, y-ray emission from unresolved star-forming galaxies, 
y-ray emission from dark matter annihilation, and y-rays from 
electromagnetic cascades initiated by very high and ultrahigh 
energy particle interactions with the extragalactic background 
light. The spectra of these potential contributions to the EGB 
differ considerably from that of the FSRQs (Silk & Srednicki 
1984; Stecker et al. 1985; Rudaz & Stecker 1988; Stecker 
& Tylka 1989a, 1989b; Rudaz & Stecker 1991; UUio et al. 
2002; Ando et al. 2007; Kalashev et al. 2009; Siegal-Gaskins 
& Pavlidou 2009; Berezinsky et al. 2011; Ahlers et al. 2010; 
Venters 2010). However, recent Fermi observations have placed 
significant constraints on dark matter annihilation (Cirelh et al. 
2010; Abdo et al. 2010a; Ackermann et al. 2010a), and presently 
there is no clear evidence of annihilation features above the 
background continuum. As such, it appears that any putative 
contribution to the EGB from dark matter annihilation is 
relatively minor. The possible contribution to the EGB from 
electromagnetic cascades is constrained by the relative steepness 
of the EGB spectrum, though cascades could play a role at 
higher energies (Kalashev et al. 2009; Berezinsky et al. 2011; 
Ahlers et al. 2010; Venters 2010). An apparent explanation for 
the featureless power-law spectrum of the EGB as presently 
deduced could be that unresolved blazars provide the dominant 
contribution to the EGB, given that their collective spectrum is 
roughly consistent with that of the EGB. 

Therefore, we conclude that, contrary to the result given by 
Abdo et al. (2010k), the Fermi observations do not rule out 
the possibility that the EGB is dominated by emission from 
unresolved blazars. 
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